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Abstract. We study theoretically the quantum dynamics of nitrogen molecules (N 2 ) 
exposed to intense and ultrafast x-rays at a wavelength of 1.1 nm (1100 eV photon 
energy) from the Linac Coherent Light Source (LCLS) free electron laser. Molecular 
rate equations are derived to describe the intertwined photoionization, decay, and 
dissociation processes occurring for N 2 . This model complements our earlier 
phenomenological approaches, the single-atom, symmetric-sharing, and fragmentation- 
matrix models of J. Chem. Phys. 136, 214310 (2012). Our rate-equations are used to 
obtain the effective pulse energy at the sample and the time scale for the dissociation 
of the metastable dication Nj ■ This leads to a very good agreement between the 
theoretically and experimentally determined ion yields and, consequently, the average 
charge states. The effective pulse energy is found to decrease with shortening pulse 
duration. This variation together with a change in the molecular fragmentation pattern 
and frustrated absorption —an effect that reduces absorption of x-rays due to (double) 
core hole formation—are the causes for the drop of the average charge state with 
shortening LCLS pulse duration discovered previously. 
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1. Introduction 

With the development of x-ray lasers, especially the state-of-the-art x-ray free electron 
lasers (FELs), the exploration of the interaction of intense and ultrafast x-rays with 
matter has become a scientific frontier ms\. The extreme pulse characteristics of 
operational x-ray FELs—such as the the Linac Coherent Light Source (LCLS) [ T91120] 
in Menlo Park, California, USA—are ultrahigh brightness, femtosecond pulse duration, 
short wavelengths, and transverse coherence. The novel properties of x-rays enable 
the study of ultrafast multiphoton interaction of x-rays with matter and spawned the 
subfield of x-ray quantum optics EU- Specifically, the intense x-rays from FELs induce 
two-photon absorption [22] (extreme ultraviolet (xuv) [23]), Rabi flopping [[241125] . 
photoelectron holography ra. single-pulse pump-probe experiments E3, stimulated 
x-ray Raman scattering [28], and four-wave mixing [28] (xuv [29]), suppress the 
branching ratio of Auger decay (50) or produce interference effects in Auger decay [311132] . 
Saturable [33] and frustrated absorption [l]0ii34] (xuv [35]) reduce radiation damage 
of the sample and may thus be beneficial for diffraction experiments at FELs [36]. The 
combination of intense and ultrafast x-rays with an optical laser offers perspectives for 
the control of x-ray lasing [37], high-order harmonic generation in the kiloelectronvolt 
regime [38l l4l] . and the production of high-energy frequency combs 12.43|. 

Initially, atoms were studied in intense and ultrafast soft x-ray pulses from LCLS [4]. 
Neon was chosen as a prototypical atom [4] and it was found that, for sufficiently 
high photon energies, all atomic electrons are ionized by one-x-ray-photon absorption 
and atoms even become transparent at high x-ray intensity due to rapid ejection of 
inner-shell electrons. Ionization by x-rays [4] (xuv [441145] ) progresses from the inside 
out in contrast to strong-field ionization with intense optical lasers where ionization 
of the outermost electrons occurs, i.e., there ionization proceeds from the outside in. 
The quantum dynamics of neon atoms in intense and ultrashort x-ray radiation was 
successfully described theoretically with a rate-equation model [UI16J foreshadowing 
extension to more complex systems. 

Understanding the response of molecules to intense and ultrafast x-rays poses 
substantial complications over the description of atoms [MET] Ell 07]; namely, in 
addition to photoionization and intraatomic decay processes, for molecules, there are 
also molecular fragmentation, and sharing of the charges in the valence shells among 
the nuclei which need to be regarded. The quantum dynamics induced in molecules by 
soft x-rays [IH3] [27l[34ll47 j from FELs (xuv [48H50] L has been studied in experiments 
during the last few years. The first ever study of a molecule, namely nitrogen (N 2 ), in 
intense and ultrafast x-rays was performed at LCLS [1] where N 2 was chosen because 
it is an important but simple molecule. The experiment used x-rays with a wavelength 
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of 1.1 nm, i.e., 1100 eV photon energy. We found [I] that the absorption of x-rays is 
frustrated, i.e., pulses with a comparable energy but decreasing duration lead to smaller 
average charge states compared with longer pulses and thus a reduced absorption of 
x-rays. Unlike the XUV photon energy regime in which multiphoton absorption causes 
single or multiple ionization [^8b50] . predominantly only one-photon ionization occurs 
in the x-ray regime [Hll., 2?', 31] which is a substantial simplification of the problem. 

To find a theoretical description of N 2 in intense and ultrafast x-rays, we devised 
a series of phenomenological models of increasing sophistication, namely, a single-atom 
model, a symmetric-sharing model, and a fragmentation-matrix model [34]. It was found 
that a single-atom model—which assumes that the molecular ion yields are the same 
as the ion yields that are obtained from a single atom—is not capable to describe the 
experimental data for N 2 due to a redistribution of valence charges over both N atoms 
in the molecule after ionization and prior breakup. Therefore, a symmetric-sharing 
model was tried next in which the molecular charge is distributed evenly between the 
two atoms. The model was found to lead to a similarly deficient description as the 
single-atom model. Thus the insights gained from these two models were used to devise 
heuristically a fragmentation-matrix model in which there is only a partial redistribution 
of charges between the two atoms. This model clearly reveals the relevance of the 
redistribution of charge and a different weighting of electronic processes and nuclear 
dynamics on short and long time scales for similar nominal pulse energy but varying 
pulse duration. 

In this article, we study sequential multiple ionization of N 2 induced by intense 
and ultrafast x-rays from LCLS with comparable nominal pulse energy but different 
pulse durations. In the theory section [2] we extend our modeling from [[34] by devising 
a molecular rate-equation formalism^ to describe the quantum dynamics induced by 
x-ray absorption. The molecular rate equations are formulated for neutral N 2 and 
molecular cations with single core holes (SCHs), double core holes (DCHs) on a single 
site (ssDCH), and double core holes on two sites (tsDCH) [[TQlf34lf5TT[53] . Further rate 
equations for valence ionized molecular configurations and even triple core holes round 
off the set. In total, we consider all possible one-photon absorption and Auger decay 
channels for molecular charge states up to N^" 1 " [34[j54j. We treat molecular charge 
states beyond N^ + as fragmented into single atoms and we take all single atom charge 
states into account in terms of atomic rate equations [311154] • By considering the finite 
lifetime of the metastable molecular dication N^" 1 " and the sharing of the valence electrons 
between the two nuclei of the molecule, the experimental ion yields and the average 
charge state are well reproduced in the results and discussion section [3] Our present 
molecular rate-equation model shows that the differences in the ion yields observed for 
pulses of varying durations are due to the competition between photoabsorption, Auger 

§ The model for N 2 in intense and ultrafast x-rays of Ref. [Tj is also based on molecular rate equations 
and was developed by Oleg Kornilov and Oliver Gessner. The rate equations in this work were devised 
independently and their results are consistent with the phenomenological models of I34j and further 
analysis of the experimental data in [27] unlike the previous attempt in Ref. [I]. 
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Figure 1 . Schematic of the interaction of N 2 with intense x-rays. The molecular 
processes displayed are photoionization by one-x-ray-absorption “<r”, Auger decay of 
core-ionized states “T”, dissociation of the metastable molecular ions with a 
fragmentation rate “7” , and instantaneous fragmentation “00” . In the last step of 
the interaction “00” occurs, if no core holes are present in N| + ; otherwise “cr” or ‘T” 
take place. 


decay processes, and molecular fragmentation. Frustrated absorption is explained as 
a reduction of photoabsorption by core hole formation. The effective pulse energy at 
the sample and the rate of the dissociation of the molecular dication are obtained with 
the molecular rate-equation model by comparing theoretical and experimental data. 
Our rate-equation model represents an improvement over the fragmentation-matrix 
model [M] as it describes the quantum dynamics of the involved molecular charge 
states and the breakup process. Nonetheless, the results from our model agree well 
with our earlier findings [34] . Conclusions are drawn in section [4] All details of the 
molecular rate-equation model and the calculations in this article are provided in the 
Supplementary Data 0- 

Atomic units [56, 57] are used throughout unless stated otherwise. For the 
conversion of a decay width Y in electronvolts to a lifetime r in femtoseconds, we use 
the relation r = f = a ™ 2eVfs . 

2. Theory 

We develop a rate-equation description for N 2 in intense and ultrafast x-rays which 
includes the molecular processes schematically depicted in figure [0 Initially, all 
N 2 molecules are in their groimd state. The absorption of an x ray creates the molecular 
cation N^" which is frequently metastable with respect to dissociation. Thus remains 
a molecule until either another x ray is absorbed or, if there is a core hole, Auger decay 
takes place [58H6T] . Both processes promote to the dication N^ + . Even that, 
in our case, is predominantly produced by Auger decay, is frequently metastable with 
a comparatively long lifetime of ~100fs with respect to dissociation [62][63]. There 
are three channels for Nj" 1 " to evolve into; first, N^" 1 " fragments causing the end of 
the molecule; second, an x ray is absorbed, and, third, if there is a core hole, Auger 
decay occurs. Both the second and the third channel promote N^ + to N^ -1 ". If N;] + 
has no core holes, we assume that it fragments immediately into atoms; otherwise, 
again x-ray absorption and Auger decay may occur turning N^ -1 " into Ng which is 
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considered to fragment immediately into two N 2+ atoms [631165] . From here onwards, 
only atomic fragments are considered and their further evolution is governed by atomic 
rate equations. 

Molecular rate equations are formulated for N 2 incorporating the steps mentioned 
in the previous paragraph. Altogether there are thirty-three molecular configurations 
from neutral N 2 to triply charged with different combinations of the electronic 
configurations of the two atomic sites. The quantum dynamics of the atomic fragments 
after breakup is described with an atomic rate-equation model which is modified in order 
to funnel probability from the molecular rate equations into the atomic rate equations. 
All in all, we extend nine of the thirty-six atomic rate equations to describe the transfer 
of probability from molecular configurations to atomic configurations. The molecular 
and atomic rate equations incorporate all energetically accessible atomic one-photon 
absorption processes together with fluorescence and Auger decay of multiply-charged 
nitrogen atoms [53] : simultaneous multiphoton absorption [22l1361166] and multielectron 
processes such as shake off ra and double Auger decay ra, however, are neglected. 
In what follows, we write down selected molecular rate equations to show the guiding 
principles. The complete set of equations is provided in the Supplementary Data 0- 

2.1. Neutral nitrogen molecule N 2 

We specify probabilities Pijk/mnif) at time £ for 0 < < 2 A 0 < k,n < 

3 to find an N 2 molecule in which the individual N atoms have the electronic 
configurations 1 s l 2s^2p k and ls £ 2s m 2p n , respectively, which are abbreviated by {zj 
and {£mn}. Molecular configurations for N 2 are written succinctly as {ijk,£mn 
X-ray absorption by N 2 is governed by a single rate equation: 

dP 2 23,223<X) _ o„ 6 T /4\ t-w 

--j^-— —2 cr 223 -r 22 3 i223 (£) J x (£) , (1) 

with the total one-x-ray-photon absorption cross section ay,*, of an N atom in 
configuration {ijk} where i = 2, j — 2, k = 3 denotes the neutral N atom in the 
ground state. The factor of 2 accounts for the two atoms in N 2 and the minus sign 
indicates a depletion of N 2 in the ground state by x-ray absorption. Further, J x (£) is the 
x-ray photon flux at time t. Note that the x-ray flux has, in principle, both temporal and 
spatial dependence due to the longitudinal and transverse profiles of the x-ray beam [34]. 
For simplicity of notation, we do not explicitly write spatial coordinates for the quantities 
involved but state all equations for a flux that depends only on time J x (£). 

2.2. Molecular nitrogen cation Nj 

For the singly-charged molecule Nj with a SCH, the governing rate equation reads 

——~7T~—~ ^ 2 P 223 ,223(0 Mt) (2) 



|| Note that {ijk,imn} = {£mn,ijk} holds, i.e., only distinct molecular configurations up to a 
permutation are considered to avoid double counting. 
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— (°'123 + cr 223) ^123,223 if) Jx(t) ~ Ti23 Pl22,,222,{t) j 

where <Je mn< ^ijk stands for the one-x-ray-photon absorption cross section which causes 
an N atom with the configuration {ijk} to transition into the configuration {£mn}. 
With Tjjfc we denote the total decay rate of an N atom in configuration {ijk} to any final 
state. The first term on the right-hand side of (J 2 j) is the rate of production of {123, 223} 
by x-ray absorption; the second term describes the rate of depletion by x-ray absorption; 
and the third term represents the rate of decay. 

Further molecular rate equations are formulated for the other cationic 
configurations of N}“ which are {213,223} and {222,223}; the resulting molecular rate 
equations are expressed analogously to o m Only few nondissociated N}~ molecules 
were detected in the experiment [I] which justifies that we do not treat them as a channel 
but break them into atomic fragments at t —» oo. 


2.3. Molecular nitrogen dication N 2+ 


For the molecular dication N| + with only valence vacancies, we have the exemplary rate 
equation 


dd^2l2,223 if'} 
(it. 


0212<-213 ^213,2230) Jx(t) + 0212-<-222 P 222,223(f) Jx(t) 

~ (0212 + 0^223) -^212,223(f) Jx(t) + h212<—123 -Pl23,223(f) 

— 7 ^212,223(f) , 


(3) 


where Te mn ^ijk is the decay rate of the atom in configuration {ijk} into {£mn}. The 
fragmentation rate of the valence-ionized molecular dication N^ + is indicated by 7 which 
is presumed not to depend on the molecular configuration. We assume that all N^ + 
eventually dissociate because only few nondissociated N^ + molecules were detected in 
the experiment pjj. 

Since N^ + dissociates either symmetrically into N + + N + or asymmetrically 
into N + N 2+ , the parameters / N + +N + and /n+n 2 + are introduced [31] to account for 
the probabilities to dissociate into the fragmentation channels 


N 2+ /n++n+ N + + N + ? ( 4a ) 

N 2 + /n+n2+ n + N 2+ ^ (4b) 


with / N +_|_ N + + /n+n 2 + = 1- The molecular probability of N^ + upon breakup needs to 
be funneled into the atomic rate equations. For example, the atomic rate equation of a 
neutral N atom has to be modified to 
dT*223 (f) 


d t 


— — 0223 -P223(f) Jx(t) 


+ /n+n 2 + 


7 


-P203,223(f) + T212,223(f) + -Pm,223(f) 
-P213,213(f) + ^213,222(f) + ^222,222(f) 


(5) 






Rate equations for nitrogen molecules 


7 


where Pe mn (t) is the probability to find an N atom at time t in electronic 
configuration {imn\. The last term on the right-hand side of ([5]) is new compared with 
the original atomic rate equation and it is responsible for funncling probability of neutral 
N atoms due to molecular fragmentation from the molecular configurations {203,223}, 
{212, 223}, {221,223}, {213, 213}, {213,222}, and {222, 222}. Similar modifications are 
made to the atomic rate equations for the configurations {213}, {222}, {203}, {212}, 
{221}, {113}, {122}, and {023} [H]. Other atomic rate equations are not altered. 

For a valence hole and a SCH on one atom, the dynamics is quantified by the 
molecular rate equation 


dPi 


113,223 


(t) 


d t 


for a tsDCH by 


dPi 


123,123 


(t) 


dt 


and for a ssDCH by 


dP r 


023,223 


(t) 


d t 


— & 113<—213 P213,223(X) ^x(t) + <Xll3«-123 Pl23,223(A Jx(t) (6) 

— (CTH3 + 1 ^ 223 )Pll3,223 if') Jx(t) ~ Tn3 Pl 13,223H) j 

— &123<—223 Pl23,223 if) Jx{t) ~ 2 Cr 12 3 Pl23,123 (A Jx{t) (7) 

— 2 r i23 Pi23,l23 if) , 

= ^023^-123 Pl23,223 if) Jx{t) ~ { a 023 + <7223^023,223 (A Jx{t) (8) 

— r023 Po' 


023,223 if') • 

We see from these expressions, that molecular fragmentation is not included when a 
SCH or DCH is present because the time scale is, in this case, determined by the decay 
time of the vacancies which is much shorter than the fragmentation time. 

The other dicationic configurations of N| + are {213,213}, {222,222}, {123,213}, 
{123,222}, {213,222}, {203,223}, {221,223}, and {122,223}; the involved molecular 
rate equations are formulated in analogy to the above scheme [5~4| . 

2-4- Molecular nitrogen trication N^" 1 " 

Further photoionization or Auger decay of the molecular dication N^" 1 " produces the 
trication N^ + . No bound states were found for N} + ions in the midst of many Coulomb 
repulsive potentials [68]. Therefore, we assume that those N^ -1 " ions with no core 
holes fragment into N 2+ and N + at the instant of their formation by Auger decay or 
photoionization. The probability for the formation of N^ + ions without core holes is 
small due to small valence ionization cross sections at the x-ray photon energy [ 69] . 
In contrast, N^" 1 " with core holes is still treated as a molecule, i.e., no immediate 
fragmentation is assumed; the N^" 1 " are only broken up by Auger decay or photoionization 
because of the faster time scale of Auger decay compared with molecular fragmentation. 

For the triply charged molecular configurations, with a SCH and two valence holes, 
we have the exemplary rate equations 
dPl21,223 (t) 


d t 


— C r 121<-221 P221,223 if) Jx(t) + 121^—122 Pl22,223(A Jx(t) 
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(0121 + ^ 223 ) ^121,223 if) Jx(t) + ri21<-023 -^023,223 (t) (9) 

r m Pi21,223(1) j 


d P, 


1(0 


— 0123-5-223 ^212,223^) Tx(t) + cr 212<-213 -^123,213^) Jx(t) (10) 
+ 0212-^222 -fl23,222 W ^x(^) — ( cr 212 + 0123) ^212,123(0 Tx(t) 


d£ 


2 ^12-5-123 -Pl23,123 if) ~ Tl23 -^212,123 


(*), 


a tsDCH and a valence hole 
dAl3,123(^) 


d£ 


— 0113-5—213 -fl23,213(^) Jx(t) + 0123^-223 -^113,223^) </x(£) 

+ 2 0n 3< _ 1 23 -^123,123(^) Tx(t) — (0113 + 0123) 
x Pi 13,123(0 Jx(t) — (r il3 + r i23 ) -Pil3,123(0 j 


( 11 ) 


a ssDCH and a valence hole 
d-Pm 3 223 ^) 


— 0013-5—113 -fll3,223(0 Jx(t) + 0013-5—023 -^023,223 if) Jx(t) 
~ (0013 + 0223) -^013,223 if) Jx(t ) — r 0 i 3 -^013,223(£) j 


d t 


( 12 ) 


and a ssDCH and a SCH 

--—- = 2 0023-5-123 Pl23,123(t) Jx{t) + 0123-5-223 -^023,223^ ) Jx(t) (13) 

at 

— (0023 + 0123) -^023,123 if) Jx(t ) — (r 0 23 + Ti 23 ) -^023,123(0 • 

The other tricationic configurations of N| + are {113,213}, {113,222}, {122,213}, 
{023,213}, {203,123}, {112,223}, {103,223}, {122,222}, {023,222}, {221,123}, 
{122,123}, and {022,223}; the resulting molecular rate equations are written 
accordingly m- 


2.5. Higher charge states of a nitrogen molecule 

Molecular effects are centered around the steps described in section 12.11 12.21 12.31 
and 12.41 [34]. For even higher charge states, the interaction is basically those of 
independent atoms and no molecular rate equations are formulated. Instead, immediate 
fragmentation is assumed for the molecular tetracation N 2 + where the four charges are 
shared equally between the two nitrogen nuclei [641165] . The probability of the fragments 
is funneled into the atomic rate equations which govern the time evolution from this 
point onwards [54] . 


3. Results and discussion 

To predict the ion yields and the average charge state of N 2 in intense and ultrafast 
x-rays, we use a number of computational parameters which were either taken from 
experiments mm or calculated [34] • The x-ray field parameters used here are listed 
in table [D and are the same as the ones in [T1I34}. The expression for the beam profile 
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of LCLS pulses can be found in [34]. The dissociative photoionization cross sections 
of a SCH in N 2 are known from experiments at third-generation synchrotrons [61] . 
From these we find [34) for the fragment ion ratios in (j4J) the values /n++n+ = 0.74 
and / N+N 2 + = 1 — / N++N+ = 0.26; we use them to parametrize our rate equations to 
ensure that the ion yields approach the values of SCH decay from synchrotrons in the 
limit of low x-ray intensities and fluences. Auger and fluorescence decay widths and 
the corresponding transition energies, electron binding energies, and one-x-ray-photon 
absorption cross sections for multiply-ionized N atoms can be found in [34] where they 
were determined with ab initio computations. The joint molecular and atomic rate 
equations form a system of ordinary first-order linear differential equations which is 
solved numerically with Mathematica [54J[55] with the initial condition that N 2 is in its 
ground state, i.e., -P 223,223 (— 00 ) = 1 , and all the other molecular states and all atomic 
states are unpopulated, i.e., Pijk,imn(—oo ) = 0 with 0 < i,j,£,m < 2 A 0 < k,n < 
3 A -i(i = j = £ = m = 2Ak = n = 3) and Pijk(-oo) = 0 with 0<i,j<2A0<k<3. 

The ion yields of atomic fragments from molecules with a charge of 1 < j < 7 are 
the molecular ion yields 


Yj = 

J 1 


Pi 


Pn 


(14) 


which are experimentally mensurable quantities and are defined as the renormalized 
charge-state probability where 


Pj ^ ^ Ptmn{oo) ; 

£+m+n=7—j 
0<£,m<2 A 0<n<3 


(15) 


is the probability to find an atomic fragment with charge 0 < j < 7 for t —> 00 when 
molecular breakup is always assumed. From the ion yields (114]) . we calculate the average 
charge state via 

7 

9 = ( 16 ) 

3 = 1 

it is an indicator of the amount of charge found on molecular fragments. 

In the experiment, the beam line transmission is not particularly well determined 
due to beam transport losses, which leads to a large uncertainty in the pulse energy 
actually reaching the sample; we estimated that only 15%-35% of the nominal pulse 
energy arrived at the sample [111271134] . Therefore, the effective pulse energy Ep and the 
rate of molecular fragmentation 7 [Equation (j3J) and ([5])] in our calculations are unknown 
and need to be found by comparing theoretical results with experimental data. For this 
purpose, we jointly determine Ep and 7 by minimizing the criterion 


C'j(-£/p,y ) |(7expt (Ztheo(A'p, 7 ) | T 16 [FTxpt j Y theo j(Ep, 7)| , (1*0 


it comprises the absolute value of the difference between the experimental g expt and 
the theoretical qtheo(Ep,A) average charge state and the absolute value of 16 times the 
difference between the experimental Y expt j and the theoretical Y theo j(E P , 7) ion yield of 
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Table 1 . Parameters of the molecular rate-equation model. The LCLS nominal 
FWHM pulse duration is tx and the effective pulse duration is [M] . A nominal pulse 
energy of 0.15 mJ is specified for 4fs pulses and 0.26 mJ for the remaining three pulse 
durations of 7 fs, 80 fs, and 280 fs. The effective LCLS pulse energy from molecular rate 
equations is Ep, i.e., the energy which arrives at the sample. The fragmentation time 
for nominal pulse durations from an optimization (ED with respect to the N + yield 
is 7 ^+ and with respect to the N 2+ yield it is 7 ^ 2 +; correspondingly 7 ^ 7 1 an< l 7^7+ are 
the fragmentation times for effective pulse durations. No fragmentation times could 
be determined for 4fs and 7 fs pulses. The LCLS photon energy is 1100 eV. 


T X [fs] 

T k [ f s] 

Ep [mJ] 

7 n + [fs] 

7N+ 1 [fs] 

7 N 2 + [fs] 

7n 2+ [fs] 

280 

112 

0.38 x 0.26 

246 

95 

201 

76 

80 

40 

0.30 x 0.26 

83 

40 

60 

26 

7 

2.8 

0.16 x 0.26 





4 

1.6 

0.23 x 0.15 






either N + for j — 1 or N 2+ for j = 2. The second summand in (fT7|) causes the criterion 
to have a fairly sharp global minimum where the prefactor of 16 is somewhat arbitrary 
but chosen large enough to ensure that the second summand makes a large contribution 
thus effecting a good agreement between experimental and theoretical average charge 
states provided that simultaneously also the chosen experimental and theoretical ion 
yield agree well; only N + or N 2+ ion yields (J4j) are considered in the second summand 
because only these are influence directly by 7 [Equation (J3]) and (j5J)]. The ion yields 
from higher charge states are only indirectly impacted by fragmentation [34] . 

We list in table |T| the effective pulse energies and fragmentation times obtained 
by minimizing the criterion (fT7D for 80 fs and 280 fs nominal LCLS pulse durations. 
It was found that the nominal pulse durations specified by the accelerator electron 
beam parameters ina are substantially longer than the effective pulse durations at 
the sample [1 ED] for which values for the fragmentation time are also given. As the 
dependence of 4 fs and 7 fs pulses on 7 is very weak, it was not possible to extract it 
from (1X71) . In order to obtain the effective pulse energy for the short pulse durations, 
7n+ is used with the value from the nominal duration of the 280 fs pulses. The 
percentages of the nominal pulse energies that effectively arrive at the sample lie almost 
in the range of 15%-35% which was specified in the experiment [T]. The percentages 
in table [Q of this work have a close agreement with those of the fragmentation-matrix 
model in Table I of [34] which are 26%, 16%, 25%, and 31% for pulse durations of 4fs, 
7 fs, 80 fs and 280 fs, respectively. This means that the effective energy of the LCLS 
pulses arriving at the sample decreases with the shortening of the pulse. A reanalysis of 
the experimental data of [l] in [27] confirmed the predicted drop in the effective pulse 
energy. 

Comparing 7 ”+ and 7 n2 ' + for a specific pulse duration in table |T] reveals a noticeable 
difference which is due to the fact that different pathways contribute to these ion 
yields. For example, N + ions are not only produced by (14X1) but, e.g., also by valence 
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ionization of N; and N 2+ ions do not only result from (14bj) but, e.g., also from Auger 
decay of a SCH in a N + ion. As neither N + nor N 2+ are solely produced by molecular 
fragmentation, this causes a systematic error, if we use (JED, which can be estimated 
from the differences between the respective fragmentation times in table |T] for the same 
pulse duration. Another source of error are the variations of the fragmentation times 
for a specific ion yield with the pulse duration, e.g., 7 ^+, obtained from 80 fs and 
280 fs pulses. This happens because the contribution of the pathways varies when the 
pulse duration is shortened as the peak x-ray intensity increases. The same holds true for 
the fragmentation times found from the effective pulse durations which are not known 
particularly precisely. In contrast to the theoretical results from the fragmentation- 
matrix model [33, which were found to depend on the pulse duration only weakly, for 
molecular rate equations the determination of the fragmentation rate 7 is sensitive to 
the actual pulse duration. Overall, we find 7 -1 to be in the range of 40-250 fs which is 
comparable to the value of ~100fs found in [63]. The fragmentation time 7 _1 is much 
longer than the lifetime of a SCH of 6.73 fs [34j but is faster than the time for SCH 
production by photoionization for the 280 fs pulse which is 348 fs at the peak intensity. 
For the 4fs short pulse, the lifetime of a SCH is comparable with the time between two 
photoionizations of 10.1 fs whereas the time of molecular fragmentation 7" 1 is about an 
order of magnitude longer [54] - 

The ion yields of N 2 from our molecular rate-equation model and the experiment 
are depicted in figure [ 2 ] where also theoretical ion yields are shown when certain 
DCH channels are closed. The results from the molecular rate equations agree nicely 
with the experimental data when all channels are included. Considering the results for 
which the DCH-containing channels are omitted, i.e., only the SCH channel is included, 
we see that the dominant features of the ion yields are reproduced. Specifically, the 
influence of DCHs is small for 280 fs pulses but DCHs have a substantial influence on 
the ion yields for the 4 fs pulse. This result is consistent with the analysis in Section V 
of [23]. The longer the pulse duration, the lower is the peak x-ray intensity of the 
pulse. Hence the rate of photoionization is fairly small for long pulses and Auger 
decay and molecular fragmentation are the dominant processes compared with DCH 
formation. On the contrary, for short pulses, the molecules are ionized faster resulting 
in DCHs but nuclear dynamics can be neglected because of the short pulse duration 
and valence charges are mostly shared equally between the two nitrogen atoms [34| . 
The artificial closing of DCH channels inhibits x-ray absorption and thus leads to an 
underestimation of the amount of ionization which means that the ion yields of low 
charge states in figure [ 2 ] are higher than the ones from the complete model and the 
higher charge states are lower, respectively. 

Merits of the molecular rate-equation model over the fragmentation matrix model 
from [34] are that the former takes explicitly the time scale of molecular fragmentation 
into account and offers a description of the dynamics of the population of charge 
states whereas the latter provides only fragmentation constants. In order to reveal 
the important role of the fragmentation rate introduced in the molecular rate-equation 
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Figure 2. (Color) Ion yields of N 2 ionized by LCLS x-rays with FWHM durations of 
(a) 4fs, (b) 7 fs, (c) 80 fs, and (d) 280 fs. Experimental data (black bars); results of the 
molecular rate-equation model with all channels included (red bars), only with SCH 
and tsDCH channels included (green bars), and only with SCH channels included 
(blue bars). The fragmentation time of the molecular dication N| + is y^ 1 = 246 fs for 
all pulses but 80 fs pulses for which y _1 = 83 fs is used. Other parameters are taken 
from table [U 


model [Equation (j3J) and (|5j)]. in figure [3], a comparison of the ion yields is shown which 
are obtained from our molecular rate-equation model in the two extreme cases of 7 —> 00 
and 7 = 0 together with the results of the single-atom model and the experimental 
data. When immediate fragmentation at the instant of the formation of the metastable 
molecular dication N^ -1 " is assumed i.e., 7 —y 00 , the molecular rate-equation model 
still gives similar ion yields as the experimental data. Specifically, the theoretical ion 
yields follow the same trend as the experimental ones. Yet we see clearly from figure [3l 
that the rather good agreement for 4fs pulses becomes progressively worse for longer 
pulses. This observation reveals the increasing influence of the fragmentation time scale 
which eventually becomes comparable with the FWHM duration of the long pulses. As 
molecular breakup occurs immediately in (J3J) and (J5j) and the hierarchy of molecular 
rate equations is quenched at ([3]) because all probability is funncled into the atomic rate 
equation (j5j), several of the molecular rate equations for the trication are thus not used. 
Conversely, the other limit of the fragmentation time scale is 7 = 0 which means that 
fragmentation of the metastable molecular dication N| + is neglected. One can see clearly 
from figure [3] the large differences between the experimental results and the theoretical 
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Figure 3. (Color) Ion yields of N 2 ionized by LCLS x-rays with FWHM durations of 
(a) 4 fs, (b) 7 fs, (c) 80 fs, and (d) 280 fs. We display experimental data (black bars) 
alongside theoretical results of the molecular rate-equation model with all channels 
for 7 —> 00 (red bars) and for 7 = 0 (green bars). Further, we give theoretical 
results of the atomic rate-equation model (blue bars). Other parameters are taken 
from table [T] 


calculations. The ion yields in this limit agree well with the atomic rate-equation model 
but poorly with the experimental data. Specifically, the molecular rate equations in 
this limit do not describe the experimental observation that the ion yields of the lowly 
charged ions are higher than those of the highly charged ions (see [Mj for details of the 
single-atom model). We also observe in figure [3] that the variation of the ion yields with 
the pulse duration is substantially reduced for the atomic rate equations which indicates 
that the heightened sensitivity to the pulse duration is due to molecular fragmentation. 
In other words, the metastability of the molecular dication N^ + plays an essential role 
in the interaction of N 2 with x-ray pulses and its finite rate of dissociation must not be 
neglected. 

I 11 figure SI we display the average charge state q of (TTBD obtained from the 
experiment [T] and the molecular rate-equation model. Figure [4] shows that the model 
reproduces very well the decreasing tendency of q with the decrease of the pulse 
duration. This comes to us as no surprise because we fitted on the theoretical q 
to the experimental value in order to determine the effective pulse energy at the sample. 
Yet the results for which certain DCH channels are closed provide novel insights. 
Inspecting figure SI we find that the SCH channel determines the main character of 
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Figure 4. (Color) Average charge state q [Equation (fThl) ] of N 2 ionized with 
LCLS pulses of varying duration. The measured and computed data points are 
shown and connected by spline interpolation. The experimental q is shown as black 
circles, the theoretical q , calculated with molecular rate equations is shown as red 
squares (SCH, tsDCH, and ssDCH channels), green diamonds (only SCH and tsDCH 
channels), blue triangles up (only SCH channels), and magenta triangles down 
(constant pulse energy, SCH, tsDCH, and ssDCH channels). The parameters are the 
same as in figure [2] 



the dependence of q on the pulse duration. It describes well the nuclear and electronic 
dynamics of N 2 for long pulses. That is because for, e.g., the 280 fs pulse, Auger 
decay is the fastest process; a core hole formed by photoionization of an inner-shell 
electron is typically refilled by Auger decay before another core electron is photoionized. 
However, the restriction to only SCH channels significantly underestimates q for short 
pulses compared with the experimental result which can only be obtained by including 
the molecular configurations with tsDCH and ssDCH. This means that there is an 
appreciable probability of ionization of a second core electron before Auger decay of 
the core hole occurs. Furthermore, one can see from figure [2] and figure [4] that tsDCHs 
make a stronger contribution than ssDCHs because of the larger cross section for the 
formation of a tsDCH—as the x-ray photoabsorption cross section of an atom with no 
core hole is larger than that with a single core hole—and the lifetime of tsDCH is longer 
compared with the lifetime of a ssDCH [34] . 

Overall, there is a drop in q for short pulses compared with long pulses in figure |4] 
This observed frustrated absorption is ascribed to three effects: first, the decrease in 
the x-ray absorption cross section due to SCH and DCH formation, i.e., the effect of 
frustrated absorption cfl second, due to a predominantly symmetric sharing of the 

This is in contrast to |.‘M| where we focused on the sharing of charges due to molecular breakup on 
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charges for short pulses upon breakup as the molecular ion remains intact during the 
interaction in contrast to long pulses with fragmentation in the course of the interaction, 
and, third, the variation in the pulse energy at the sample which drops noticeably with 
shortening pulse duration [Table [TJ. We also show q in figure Q] when the pulse energy 
is constant. In this case, the drop of q is only caused by the change of the breakup 
pattern and the effect of frustrated absorption. Considering only 4 fs and 7 fs pulses, 
there ought to be no significant change of the breakup pattern as the time scale for 
molecular dissociation of the metastable N^ + is much longer than the pulse durations. 
Nonetheless, there is a pronounced drop of q for 4 fs pulse compared with 7 fs pulses 
which we ascribe entirely to frustrated absorption. 

Figure [4] of this work corresponds to figure 4 of [34] which shows the experimental q 
together with q from the single-atom model that represents an upper limit to g, the 
symmetric-sharing model which provides a lower limit to g, and the fragmentation- 
matrix model that reproduces the experimental q very well for short pulses but is less 
accurate for longer pulses as the fragmentation time scale is not explicitly included in 
the model. The agreement between the curves from the fragmentation-matrix model 
and the molecular rate-equation model is not a coincidence; it shows that the main 
physical effects are described. 

Good agreement between theoretical calculation and experimental measurement of 
the ion yields and the average charge state of N 2 for LCLS pulses with varying durations 
was also obtained in jT] saying that molecular rate equations were used. However, there, 
the molecular valence charge dynamics is modeled differently to here because there 
N 2 ” 1 " ions were assumed to fragment into N 3+ + N + or N 2+ + N 2+ with different ratios 
of N 3+ +N + versus N 2+ +N 2+ —depending on whether DCH decay or two photoionization 
with Auger decay cycles occurred—which were claimed to be adjusted to obtain the best 
agreement with the experimental data in the calculation there. Furthermore, there the 
effective pulse energy at the sample was adjusted in the range 17%-21% of the nominal 
pulse energy—which was always [table [Q taken to be 0.26 mJ in [1]—for pulses with 
FWHM durations of 4 fs, 7 fs, 80 fs, and 280 fs. Our computations reveal a stronger 
dependence of the average charge state on the pulse energy in the range 16%-38% of 
the nominal value for the pulse energy of different durations where the nominal pulse 
energy for 4fs pulses was 0.15mJ and 0.26mJ otherwise [tabled]. The differences 
between the model discussed here and the model of {l] cannot be explained by the fact 
that SASE-type pulses were used in P instead of pulses with a Gaussian temporal shape 
employed here. Namely, the analysis in [34] reveals that the impact of the spikiness 
of SASE pulses compared with a Gaussian pulse is small and may be neglected in 
good approximation. We believe that three aspects of our model cause the observed 
differences: first, all possible one-photon absorption and decay channels of an N atom 
are determined from ab initio computations [34] and are included in our model, second, 
the molecular configurations up to Ng 3 " are treated, and, third, the crucial dynamics 
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of dissociation of the metastable molecular dication N^ + is considered by introducing 
phenomenologically the fragmentation rate 7 . 

4. Conclusion 

We describe the interaction of intense and ultrafast x-rays with N 2 . For this purpose, 
we lay the theoretical foundation for the description of multi-x-ray-photon absorption 
by molecules in rate-equation approximation and investigate theoretically the quantum 
dynamics of N 2 exposed to x-ray pulses from LCLS. The results are determined by 
the competition between photoabsorption, decay processes, and molecular dissociation 
for x-ray pulses of varying duration and similar nominal pulse energy. Molecular 
configurations up to trications are included in the model taking into account the 
fragmentation dynamics of the metastable molecular dications and the redistribution 
of the valence electrons between the two atomic nuclei upon breakup. Atomic 
fragments after dissociation of the molecule are treated by atomic rate equations, 
i.e., no molecular effects are regarded anymore. The molecular rate-equation model 
describes the ionization dynamics observed and allows us to calculate theoretical ion 
yields and the average charge state which are then compared with the corresponding 
experimental quantities. Thereby, the effective pulse energy at the sample and the rate 
of fragmentation of the molecular dication are determined from a comparison of 
theoretical with experimental data. We find a substantial progression in the effective 
pulse energy which decreases pronouncedly from long to short pulses. The ion yields 
and effective pulse energies at the sample from the molecular rate equations agree well 
with those obtained previously with the fragmentation-matrix model of [34]. 

I 11 [I] frustrated absorption was observed as a drop of the average charge state 
for short pulses compared with long pulses. There the most pronounced impact on 
the average charge state is due to the variation of the effective pulse energy as our 
previous [34] and present analysis reveals. Furthermore, for short pulses the molecular 
ions remain intact and fragment mostly in terms of symmetric sharing of charges 
whereas for long pulses fragmentation occurs frequently prior further absorption of 
x-rays leading to higher charge states in the latter case compared with the former. 
The effect of frustrated absorption describes the situation that core-hole formation 
reduces the probability for further x-ray absorption. It makes a significant contribution 
for short pulses compared with long pulses as DCH channels are more important 
in the former. Specifically, we find for long pulses with moderate x-ray intensity 
that photoabsorption is the slowest process; it is slower than the dissociation of the 
metastable dications Molecular configurations with SCHs dominate where Auger 
decay quickly refills SCHs and thus increase the chance for further photoabsorption in 
the course of the interaction with the pulse because the photoionization cross section of a 
filled core shell is larger than the one of an only partially filled core shell. Conversely, for 
short pulses with high x-ray intensities, the rates of Auger decay and photoionization are 
comparable. In this case, there is a much larger probability for the production of DCHs 
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by sequential absorption of two x-rays within the lifetime of the core hole. Therefore, 
further x-ray absorption is reduced—compared with longer pulses of the same effective 
pulse energy—causing a smaller population of higher charge states. Thus frustrated 
absorption results in a pronounced drop of the average charge state in addition to the 
variation of the effective pulse energy and the variation of fragmentation patterns. 

We would like to point out that our molecular rate-equation model works well 
for x-ray pulses at a photon energy of 1100 eV. We anticipate that the model remains 
applicable for lower photon energies up to the point when the photon energy still exceeds 
the highest ionization potential in the model of 667.05 eV for the ionization of the last 
A'-shell electron of an otherwise electron-bare nucleus [34] as then resonances still play 
no significant role. Yet more experimental data are needed in order to assess whether 
our model stays valid for photon energies higher than HOOeV because then the single¬ 
electron response model may become insufficient and multielectron effects such as shake 
off processes are noticeable [4]. Similar models to the ones discussed in |34j and here 
should be used to investigate other diatomic molecules in order to find out how general 
our approach is and whether also heteronuclear molecules are well described. Also the 
role of interatomic electronic decay nn ought to be investigated in the future. We are 
confident that our theory can even be extended to molecules with more than two atoms. 
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